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0\ ■ ABSTRACT 

^: 

^2 I Temporal analyses of the prompt gamma-ray and X-ray light curves of gamma-ray 

^ I bursts reveal a tendency for the burst pulse time scales to increase with decreasing 

CN ' energy. For an ensemble of BATSE bursts, Fenimore et al. (1995) show that the energy 

dependence of burst peak durations can be represented by At oc E~'^ with 7 ~ 0.4- 
0.45. This power-law dependence has led to the suggestion that this effect is due to 
radiative processes, most notably synchrotron cooling of the non-thermal particles which 
C3 ■ produce the radiation. Here we show that a similar power-law dependence occurs, under 

I certain assumptions, in the context of the blast-wave model and is a consequence of 

I the deceleration of the blast-wave. This effect will obtain whether or not synchrotron 

' cooling is important, but different degrees of cooling will cause variations in the energy 

^ ' dependence of the peak durations. 

o : 

|H , Subject headings: gamma rays: bursts — radiation mechanisms: non-thermal 

> ' 

V~j ■ 1. Introduction 

5^ , Since their discovery in the late 1960s, gamma-ray bursts (GRBs) have remained enigmatic 

despite the fact that thousands of bursts have been detected by various instruments over the last 
~ 30 years. The recent X-ray, optical and radio afterglow observations of bursts detected by 
the BeppoSAX satellite have enabled significant advances in our understanding of these objects. 
Redshift measurements associated with the afterglows of bursts GRB 970508 (Metzger et al. 1997), 
GRB 971214 (Kulkarni et al. 1998), and GRB 980425 (Tinney et al. 1998) provide convincing 
evidence that GRBs are extragalactic and may be as distant as z = 3.4 or as nearby as 2; = 0.0085. 
The temporal decay of the afterglow emission is roughly consistent with the simplest fireball/blast- 
wave interpretations (Wijers, Rees &: Meszaros 1997; Waxman 1997), and these models provide a 
specific theoretical context in which to understand the wealth of burst data which existed prior to 
BeppoSAX and which is still largely unexplained. 
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It has been argued that the spectral shapes of a substantial number of bursts, particularly 
those which have detailed spectra from the Burst and Transient Source Experiment (BATSE) 
and the other Compton Observatory instruments, are due to synchrotron emission from a shock 
accelerated distribution of non-thermal electrons (e.g., Tavani 1996). Furthermore, the spectra 
of X-ray, optical and radio afterglow emission also appear to be due to synchrotron radiation 
and exhibit characteristic signatures such as power-law behavior consistent with optically thin 
synchrotron emission in the optical and X-ray bands (Djorgovski et al. 1997; Prontera et al. 1998; 
Galama et al. 1998), self-absorption in the radio band (Katz & Piran 1997; Frail et al. 1997), and 
spectral index changes of Aa = 0.5 in the optical, indicative of synchrotron cooling (Galama et al. 
1998). 

However, even in the context of a specific dynamical and emission model, the varied complex- 
ity of the spectral and temporal properties of prompt GRB light curves remains an outstanding 
problem. In this work, we address one aspect of this problem: the tendency, for a given GRB, for 
longer burst peak durations at lower observed energy bands. We examine this issue in terms of the 
blast wave model and demonstrate how a similar tendency arises due to blast wave deceleration 
and how different degrees of synchrotron cooling can modify the basic effect to produce a range 
of behavior. In the remainder of this paper, we give a brief summary of the relevant aspects of 
the observed energy dependence of burst light curves (§ 2), describe the dynamics and emission 
properties of the basic blast- wave model and how they relate to this effect (§3), and finally, discuss 
some of the strengths and weaknesses of this interpretation and suggest some avenues for further 
investigation (§ 4). 



2. Energy Dependence of Burst Light Curves 

The general behavior of burst light curve shapes was noted by Fishman et al. (1992). The 
individual pulses tend to be sharper at higher energies and arc often characterized by fast rise times 
and slow decays, the so-called FREDs (fast rise, exponential decay). Link, Epstein & Priedhorsky 
(1993) put these descriptions on a quantitative basis by applying autocorrelation functions and 
skewness analyses to individual burst light curves. The average properties of burst peak durations 
as a function of energy were studied by Fenimore et al. (1995) for a collection of 45 bright BATSE 
bursts. They found that the average peak width can be described by At oc E~^''^^ . Furthermore, 
they found that the energy-dependent autocorrelation function possesses a nearly universal shape 
with the width of the autocorrelation function having this energy scaling. 

From studying Ginga bursts, Fenimore (1998) also notes that there is a tendency for precursor 
emission or "pre-activity" in the lower energy X-ray bands prior to the main peak at higher energies. 
He suggests that the physical processes responsible for the longer burst durations at lower energies 
also account for this pre-activity. Furthermore, Fenimore finds that the location of the light curve 
peaks are not substantially delayed at lower energies versus high. Fenimore claims that these 
phenomena appear to be inconsistent with models which only consider synchrotron cooling of 
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injected electrons since they predict longer time scales at lower energies only for the decaying part 
of the light curve and also that the light curve peaks at the lowest X-ray energies should occur only 
after the high energy emission has died away (e.g., Kazanas, Titarchuk & Hua 1998). However, 
Dermer (1998) shows that if the electrons are injected at a constant rate over a finite period of 
time, then over that time period, the pulse widths will indeed vary with energy due to synchrotron 
cooling, but the initial rise of each of the pulses will essentially line up. In this paper, we consider 
the somewhat more complex situation in which the electron injection is governed by the deceleration 
of the blast wave and the conversion of its bulk kinetic energy to non-thermal particle energy and 
magnetic fields. 



3. Blast- Wave Deceleration and the Synchrotron Spectrum 

Several descriptions of the physics of the blast wave model already exist in the literature 
and we refer the reader to Panaitescu & Meszaros (1998) and Chiang & Dermer (1998) and the 
references therein for further details. The model consists of an expanding shell of material with a 
large initial bulk Lorentz factor, F ~ 10^-10^. This shell decelerates as it interacts with an external 
medium. A shock forms, and in the shocked material, magnetic fields are generated, and particles 
are accelerated. An initial phase of nearly free expansion ends when the expanding blast wave shell 
begins to sweep up ambient matter at a rate sufficient to decelerate it significantly. During this 
deceleration phase, the bulk Lorentz factor obeys F oc where r is the radius of the shell with 
^ = 3/2 for non-radiative (or adiabatic) expansion and C — 3 for radiative expansion. 

The rate at which the bulk kinetic energy of the blast wave shell is converted to internal energy 
is given by 

— = c^A{r)mpn,^t{r)(3T{T - 1) (1) 

where A{r) is the area of the blast wave shell as a function of radius and equals 47rr^ for a spherically 
expanding shell, next is the ambient particle number density, and /3 = (1 — l/F^)^/^ (Blandford 
&; McKee 1976). In order to perform concrete calculations, it is conventionally assumed that the 
magnetic field energy density and non-thermal electron energies can be described by parameters 
and which measure the degree to which these components are in equipartition with the swept-up 
protons. In the co-moving frame of the shell, these protons have energies TmpC^. 

At the early stages of the deceleration, when the bulk Lorentz factor is F ;^ 10^, most of the 
emission seen by an observer comes from the nearest part of the blast wave shell within a narrow 
cone with opening angle ~ 1/F. Hence, the emission properties during these stages are reasonably 
well accounted for by the state of the material along the line-of-sight connecting the observer to the 
center of the explosion. We model the injected non-thermal electrons as a power-law distribution 
with a sharp lower-energy cut-off, dN/d'y oc for 7 > 7inin) all radiating from within a region of 
uniform magnetic field energy density. Fits to burst spectra indicate that the power-law index is 
typically very steep with p ^ 5 (Tavani 1996). The minimum electron Lorentz factor (and hence 
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the characteristic electron energy) is related to the equipartition parameter by 7min = £,e{fnp/fne)^ ■ 
It follows then that the peak energy of the vLjj synchrotron spectrum of this injected electron 
distribution is given by (Katz &; Piran 1997; Chiang & Dermer 1998) 

i^pcak ~ 3 X lO-^meC^n^'/t'eB^'dr^W (2) 

oc tobs-^'^/(^^+^\ (3) 

where tobs is the elapsed time as seen by the observer, and we have assumed that the burst is nearby 
(i.e., z ~ 0). We also have used the relation iobs oc r^^+^ which holds during the deceleration phase 
of the blast wave when T oc . For the discussion which follows, we will further assume that ^s, 
^e, and next are constants unless otherwise indicated. 

If synchrotron cooling and other energy loss mechanisms are not important, then the instan- 
taneous vLy synchrotron spectrum near -Epeak during the prompt burst phase can be described 
approximately by a broken power-law (Chiang &: Dermer 1998): 

vL, oc r2r4(r)(i?/i?peak)^± (4) 

oc tobs('+''^('±-'))/('^+'\ (5) 

where A+ = (3 -p)/2 for E > ^peak, A_ = 4/3 for E < ^peak (Rybicki & Lightman 1979; Katz 
1994), and in eq. |5|, we explicitly write out the iobs-dependence of uLi, at a given energy. We 
note that this time dependence only holds during the deceleration phase of the blast wave and 
only for energies less than the initial value of -Epeak (= -E'peak(ro) where Fq is the initial blast 
wave bulk Lorentz factor). The remarkably sharp peaks seen in the uFi^ spectra of a number 
of bursts are consistent with this spectrum (Tavani 1996; Schaefer et al. 1998). If synchrotron 
cooling is important, then the resulting spectrum will be somewhat more complex. In this case, 
A-|_ = (2 — p)/2, describing the softer spectrum of a cooled power-law electron distribution, and 
the i^Li, peak will be much broader and flatter, and at energies below the peak, the index will be 
A_ ~ 1/2 rolling over to 4/3 at still lower energies (Sari, Piran & Narayan 1998; Chiang & Dermer 
1998). 

Tavani (1996) and Schaefer et al. (1998) have examined the spectra of a number of individual 
bursts and find them to be consistent with an uncooled distribution of electrons (i.e., A_ = 4/3). 
However, Schaefer et al. also analyze the composite burst spectrum constructed from a sample of 19 
bright BATSE bursts, including those bursts which they also examined individually. They find an 
"average" lower energy spectral index of A_ = 0.77. Although some of the features of the composite 
spectrum can be attributed to the range of values for the peak energy among the different bursts, 
it seems likely that this softer index is due to some of the co-added burst spectra having been 
significantly affected by synchrotron cooling. 

In order to demonstrate the effect of blast wave deceleration, we consider the case of uncooled 
electron distributions. Using values C = 3/2 and p = 6, we have from eq. |5| 

UL^ oc tobs, E < ^peak (6) 

^ J. -13/4 p \ p 

OC '•obs ) ^ -C'peak) 
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and we see that the dechning portion of a blast wave pulse will only occur after -Epeak has passed 
through the energy band at which the observations are taking place. Thus, there will be an energy- 
dependent delay for the light curve decline given by inverting the expression for the energy of the 
synchrotron peak (eq. ^): 

Us oc i?-(2C+i)/4C (7) 
oc E-^-^\ C = 3/2 
oc C = 3. 

Although this analysis predicts that the timing of the pulse peaks at lower energies will also be 
displaced, in practice the intrinsic curvature of the synchrotron spectrum at i^'peak will ameliorate 
this discrepancy (cf. Fig. 1). 

Of course, synchrotron cooling will always be present in this model, though perhaps to varying 
degrees. Using our blast wave code (Chiang & Dermer 1998), we now present detailed calculations 
for cases of relatively weak (^^ = 10~^, Fq = 300) and strong {^b = 10^^, Fq = 100) magnetic 
fields. The initial bulk Lorentz factor, Fq, has been adjusted in each case so that the initial value of 
-^peak is approximately the same, and the fraction of the total kinetic energy (eq. |l]) made available 
in the form of non-thermal electrons for both cases is 10%. This latter condition ensures that 
the bulk motion of the blast wave shell in both calculations is well described by non-radiative 
deceleration (C ~ 3/2) and that any relevant differences will be solely due to the effects of the 
different magnetic fields. 

Figure 1 shows the results for the weak magnetic field case. In the upper panel, we plot the 
uLi, spectra at times tobs = lO''^^^'^'''^''''"'^ seconds after the initial fireball event. The vertical 
dotted lines indicate energies 25,57,115,320 keV — the lower bounds of the four standard BATSE 
LAD channels (cf. Fenimore et al. 1995). In the lower panel, the light curves at these four energies 
are shown, each normalized to unity at their respective maxima. We see that this simple model 
reproduces a single FRED-like pulse quite naturally. In addition, the light curve decline is delayed 
at lower energies and the energy dependence of the half-maximum pulse widths is At oc E~^'^^ 
consistent with eq. |7| for a non-radiative blast wave with no cooling. The lower energy peaks are 
also delayed, but not excessively so. We note, however, that the pulses are less sharply peaked at 
these lower energies. In Figure 2, we show the results for a strong magnetic field. As we noted 
above, the spectral peaks are much flatter and broader, and the lower energy index is close to 
A_ = 1/2. The effect on the light curves is also substantial. The light curve peaks are broader, 
but in addition, the energy dependence of the half-maximum pulse widths is significantly weaker 
with At oc E^^ '^'^. These shorter delays for the onset of the light curve declines are due to the 
combination of the blast wave deceleration and synchrotron cooling both pushing the vLy peak 
down in energy faster than either process would do so individually. 
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4. Discussion 

Given our assumptions, the effect of pulse broadening at lower energies will obtain for any blast 
wave deceleration model, whether the prompt burst emission is due to internal shocks or external 
shocks, as we have described here. Furthermore, the energy dependence of the pulse widths will 
vary depending on the degree of synchrotron cooling, and wc find that the range of this energy 
dependence (E^^'^-E^^'^^) is suggestively close to the results found by Fenimore et al. (1995). 
This simple picture, however, cannot explain the pre-activity or precursor behavior described by 
Fenimore (1998). In addition, especially for strong cooling, the pulse shapes will be less sharply 
peaked at lower energies, which implies that the autocorrelation function may not have the universal 
shape described by Fenimore et al. (1995). 

A substantially more complex blast wave model than the one we have presented here is certainly 
required to describe GRBs, and any additional complexities may either mitigate or worsen the above 
discrepancies with the observations. At a minimum, the electron and magnetic field equipartition 
parameters are not constant throughout the blast wave evolution as we have assumed. As the blast 
wave evolves, the expansion time scale also grows with increasing radius. There would then be more 
time for equipartition to obtain, and one might expect the electron and magnetic field components 
to come closer to true equipartition with the swept-up protons at later times. Such an effect may 
account for the softer pre-activity phase since -Epeak will increase as and increase until the 
temporal dependence implied by the deceleration and synchrotron losses begin to dominate and 
the effects wc have described will then control the light curve evolution. The effects of evolving 
equipartition parameters in more complex blast wave model calculations can certainly be studied, 
but one may suspect that only rather contrived electron energy and magnetic field evolution will 
reproduce the behavior described by Fenimore (1998). In this case, the ultimate answers may 
lie in doing the microphysics: magnetic field generation and particle acceleration, in other words, 
turbulent relativistic magneto-hydrodynamics. 

On the observational side, at least three additional investigations should be conducted in order 
to help guide the theoretical modeling. First, a cross-correlation or similar analysis of the prompt 
burst light curves in the various energy bands should be performed in order to quantify the energy- 
dependent delay of the pulse peaks. Second, in addition to computing properties averaged over an 
ensemble of bursts, the distributions of the energy dependences of the pulse widths and peak delays 
should also be computed. This information will determine the range of parameter space for which 
any given model must account. Finally, insofar as signal-to-noise limitations prevent time-resolved 
spectra, the spectral properties integrated over the prompt burst phase should be analyzed for 
each individual burst. Specifically, if there is evidence for synchrotron cooling in the spectrum, 
which will be signified by a z^^/^-dependence in the energy bands of interest below the uF^, peak, 
then synchrotron cooling may well be important, and in the context of the blast wave model, such 
spectral features should be correlated with the energy dependence of the pulse durations. 
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Fig. 1. — Spectral and light curve calculations in the blast wave model for weak magnetic fields. The 
initial bulk Lorcntz factor of the blast wave shell is Fq = 300, and the magnetic field equipartition 
parameter is = 10~^; other model parameters are = 1; J^ext = 10^ cm~^, p = 6. Upper panel: 
Spectra for different observer times t = IQi^i'^O-SiO. -} g. The vertical lines denote the lower energy 
bounds for the four BATSE LAD channels: 25 (solid), 57 (dotted), 115 (dashed), and 320 keV 
(dot-dashed). Lower panel: Light curves for each of the above energies. The energy dependence of 
the pulse width at half- maximum is At oc E~^-^^ . 
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Fig. 2. — Strong magnetic field case. Here = 10~^ and Fq = 100; the other parameters are the 
same as for Fig. 1. Upper panel: Spectra for different observer times t = iQi^^'^^-^fi'---} g. Lower 
panel: Light curves for the four BATSE LAD lower energy bounds. The energy dependence of the 
pulse width at half- maximum is At oc E~^'^^. 



